A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

11.09.1 NGC253; 11.09.4; 11.11.1; 11.19.2 



ASTRONOMY 

AND 
ASTROPHYSICS 



Large-scale shocks in the starburst galaxy NGC 253 



o 
o 
o 

(N 



Interferometer mapping of a ~600pc SiO/H^^CO"'" circumnuclear disk 

S.Garci'a-Burillo\ J.Martm-Pintado\ A.Puente\ R.Neri^ 

^ Observatorio Astronomico Nacional (IGN), Apartado 1143, E-28800 Alcala de Henares, Madrid, SPAIN (burillo@oan.es, 
martin@oan.es, fuente@oan.es) 

^ IRAM-Institut de Radio Astronomie Millimetrique, 300 Rue de la Piscine, 38406-St.Mt.d'Heres, FRANCE (neri@iram.fr) 
Received , accepted 



> 



o : 
o . 

o ■ 

■ 

Oh; 

6 



Abstract. This paper presents the first high-resolution 
SiO map made in an external galaxy. The nucleus of the 
nearby barred spiral NGC 253 has been observed simulta- 
neously in the v=0, J=2-l line of SiO and in the J=l-0 
line of Hi3C0+ with the IRAM interferometer, with a 
resolution of 7.5" x 2.6". 

Emission from SiO and H^^CO^ is extended in the nu- 
cleus of NGC 253. The bulk of the SiO/Hi3CO+ emission 
arises from a (600pcx250pc) circumnuclear disk (CND) 
with a double ringed structure. The inner ring, of ra- 
dius r'^60pc(4"), viewed edge-on along PA=:51°, hosts 
the nuclear starburst; the outer pseudo-ring opens out 
as a spiral- like arc up to r~300pc(20"). The kinematics 
of the gaseous disk, characterized by strong non-circular 
motions, is interpreted in terms of the resonant response 
of the gas to the barred potential. The inner ring would 
correspond to the inner Inner Lindblad Resonance (ilLR), 
whereas the outer region is linked to the onset of a trail- 
ing spiral wave across the outer ILR(oILR) . This scenario 
accounts for the unlike morphology of the maps observed 
in different molecules. 

Most notably, we report the detection of a molecular 
gas counterpart of the giant outflow of hot gas, previously 
seen in X-ray and optical lines, and tentatively identified 
as a dust chimney in the the 450/i continuum band. Two 
gas filaments appear in the SiO map to come out of the 
plane in the NGC 253 nucleus at r±60pc, delimiting the 
working surfaces of gas entrained by the hot wind. 

The SiO shows a high average fractional abundance in 
the CND of <X(SiO)>~1.5xlO-i°. This is more than an 
order of magnitude above the predicted value of a PDR. 
Moreover, X(SiO) varies at least by an order of magni- 
tude between the inner starburst region, which dominates 
the global emission, where we derive X(SiO)~l~2x 10~^°, 
and the outer region, where X(SiO) reaches a few 10~^. 
SiO abundance is also significantly enhanced in the out- 
flow (X(SiO)~3-5x lO"^*'). Different mechanisms are ex- 
plored to explain the unlike chemical processing of molec- 
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ular gas within the nucleus. Large-scale shocks induced 
by the crowding of clouds orbits across the oILR in the 
outer region and by the outflow of hot gas seem to be by 
large the most efficient mechanisms in rising the Silicon 
abundance. 



Key words: galaxies: individual: NGC 253-galaxies: ISM- 
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1. Introduction 

Among all tracers of the high density gas in spirals, the 
SiO molecule presents the most peculiar, and so far, poorly 
understood chemistry. Galactic surveys of SiO clouds show 
that the relative abundance of this species varies by 
more than five orders of magnitude (the relative abun- 
dance, X(SiO), ranges from <10~^^ to 10"''). SiO emis- 
sion is practically absent in quiescent dark clouds, indi- 
cating a high degree of depletion in grains as X(SiO)^a 
few 10~^^ (Ziurys et al 1989). In contrast, SiO is en- 
hanced in the dense and hot gas of GMCs forming massive 
stars, in particular, it appears in association with bipo- 
lar outflows where X(SiO)'-10-'^-10"^ (Martfn-Pintado 
et al 1992, Schilke et al 1997, Gueth et al 1998). These 
high abundances have been interpreted as a signature of 
shock chemistry: SiO comes into the gas phase by the 
sputtering of grains by fast shocks. At significantly lower 
abundances (X(SiO)-10-"-10-i°), SiO is predicted in 
Photon-Dominated- Regions (PDRs) , complexes of molec- 
ular gas heated by UV photons produced by OB associa- 
tions (Janssen et al 1995, Walmsley et al 1999). Finally, 
SiO emission has been reported in molecular clouds of the 
Galactic Center with X(SiO)-10"^-10~^ (Martfn-Pintado 
et al 1997, Hiittemeister et al 1998). However these clouds 
show little evidence of recent star formation. The emis- 
sion in these clouds might be connected with large-scale 
shocks rather than with star formation. Shocks are possi- 
bly driven by the bar potential or, more locally, by super- 
bubbles and supernova remnants. 
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SiO emission was first detected in NGC 253 by Maucrs- 
berger and Henkel 1991. A later study by Sage and Zi- 
urys (1995) reported the detection of J=2-l SiO emission 
in several starburst galaxies, including NGC 253. The au- 
thors concluded that the ratio I(SiO)/I(N2H-|-) is an in- 
dicator of the SFR per unit mass. However, their con- 
clusions on the origin of SiO emission are based on low 
resolution observations of only one transition. With the 
aim of studying the physical and the chemical parameters 
of SiO clouds, we mapped in the J=5-4, 3-2 and 2-1 tran- 
sitions of SiO the nuclei of three prototypical starbursts: 
M82, NGC 253 and IC342 (Martm-Pintado et al 1999, 
hereafter called MP99) . First results of our survey, made 
with the 30m telescope, underlined the different influence 
of star formation in these galaxies. First, the inferred SiO 
abundances were seen to vary from galaxy to galaxy (from 
X(SiO)~10*s in IC342 to ~10-i° in M82), showing no 
correlation with changes on the SFR per unit mass. Phys- 
ical conditions of SiO clouds, derived from the transition 
ratios, are also markedly different in our sample. 

However, any further refinement in the interpretation 
of SiO emission required a better angular resolution. In 
particular, it is crucial to elucidate which is the exact 
spatial extent of SiO emission in the nuclei of starburst 
galaxies. So far it is unclear if SiO emission is exclusively 
linked with the formation of massive stars, or alternatively, 
with the occurrence of large-scale shocks, as observations 
of clouds in our Galaxy seem to indicate. We want to 
address the question of which is the chemical scenario ex- 
plaining the abundance of SiO in galaxies. 

Fig. 1. (a,top): Emission contours of the 3mm (86.7GHz) 
continuum source at the center of NGC 253. x and y are 
offsets (in arcsec) with respect to the dynamical center 
derived in this work, at aj2ooo=00^47'"33.18*, Sj2ooo—- 
25°17'17.2"; x and y axes are parallel to the major and 
minor axes of the stellar bar, respectively {x runs parallel 
to PAhar=68°). Contours are -1.5, 1.5, 2.5, 4, 5.5, 8 and 13 
to 64mJybeam~^ by steps of 9mJybeam~^. la- noise level 
is 0.44mJybeam-^ (b,middle): SiO(v=0,J=2-l) inte- 
grated intensity contours towards the center of NGC 253. 
Contours are -0.2, 0.2 to 2.2Jy.kms~^beam~^ by steps of 
0.25Jy.kms~^beam~^. Icr-noise level in the integrated in- 
tensity is 0.07Jy.kms~^bcam~"'^. Orientation as in Figure 
la. (c, bottom): same as (b) but for the J=l-0 line of 
H^^CO+, with same contours, noise level and orientation. 

In this paper we study at high-resolution the 
SiO(v=0,J=2-l) emission in the nucleus of NGC 253, a 
highly inclined barred Sc spiral representing the archetype 
of a nuclear starburst. The bulk of its large infrared lu- 
minosity hjn ^I.GxIO^'^Lq originates in the inner 300 pc 
(Telesco and Harper 1980). The observed infrared spec- 
trum of the nucleus (Telesco and Harper 1980: Engcl- 
bracht et al 1998) as well as the radio continuum obser- 



vations made at centimeter wavelengths (Antonucci and 
Ulvestad 1988, Ulvestad and Antonucci 1991) support 
the starburst hypothesis. There is also evidence of a gi- 
ant outflow of hot gas powered by the starburst, that 
comes out of the galaxy plane (Fabbiano and Trinchieri 
1984, McCarthy et al 1987, Schulz and Wegner 1992). The 
high molecular gas mass of the NGC 253's nucleus (^2 
3xlO^M0; Canzian et al 1988; Mauersberger et al 1996) 
explains the high SFR. Water maser emission detected 
by Ho et al (1987) reveals strong star formation activity. 
The gas reservoir might have been driven inward by the 
barred potential first identified by Scoville et al 1985 in a 
K-band image of the galaxy. Several interferometer maps 
showing the molecular gas distribution of NGC 253's nu- 
cleus have been published (CO: Canzian et al 1988; HCN: 
PagHone et al 1995; CN: Hiittemeister and Aalto 1998, 
HCO+: Carlstrom et al 1990; OH: Turner 1985; H2CO: 
Baan et al 1997; CS: Peng et al 1996, hereafter P96) 
However, no high-resolution image of the dense gas distri- 
bution (n(H2)>10^cm~^) has been obtained so far, using 
an optically thin tracer. To fill this gap, we have made 
simultaneous observations in the J=l-0 line of H-'^^CO"'", 
a tracer of dense gas with low opacity. Furthermore all 
the results for SiO and H^^CO+ are discussed in paral- 
lel throughout the paper. In particular we study how the 
ratio of integrated intensities I(SiO)/I(H^^CO+) changes 
across the nucleus and discuss how this result is used to 
get the variation of the absolute abundance of SiO. 

2. Observations 

Observations of NGC 253 were made with the IRAM ar- 
ray at Plateau de Bure (France) between June and August 
1998. We simultaneously observed the (v=0,. 1=2-1) line 
of SiO (at 86.847 GHz) and the (J=l-0) hne of Hi3co+ 
(at 86.754 GHz) using the standard BC set of 4-antenna 
configurations. The 55" primary beam field of the inter- 
ferometer was phase-centered at aj20oo= 00''47™33.0^ and 
<5j20oo=-25°17'18.4". 

The SIS receivers were tuned for single side band op- 
eration at 86.731 GHz yielding system temperatures of 
100 — 200 K. The spectral correlator was adjusted to give 
a contiguous bandwidth of 420 MHz. This is equivalent to 
a velocity range of 1500 kms~^; the correlator was cen- 
tered at 86.731 GHz in order to cover both transitions. 
The effective frequency resolution was set to 2.5 MHz, or 
equivalently 8.64 kms~^ at this frequency. 

Visibilities were obtained using on-source integration 
times of 20 minutes interspersed with 4 minutes of in- 
strumental amplitude and phase calibrations. These were 
made by observing 0135-247 and 2345-167. The atmo- 
spheric phase noise on the most extended baselines ranged 
from 20° to 40°, consistent with a seeing of 0.8-1.2", typi- 
cal for summer weather conditions. The absolute flux den- 
sity scale was established on the basis of cross-correlations 
on the radio star MWC349 (with Ss7ghz = 950 mJy). 
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This is in full agreement with the measured interferome- 
ter efficiency and should be accurate to 10%. The receiver 
passband shape was determined on 3C454.3 and its accu- 
racy is better than 5% throughout the observing run. 

Cleaned maps were obtained from the visibilities using 
the standard IRAM package. The maps were 256 x 256 
pixels in extent, with a pixel size of 0.4". The synthesized 
beam, as determined by fitting a Gaussian to the dirty 
beam, is 7.5" x 2.6", oriented north-south (PA=18G°). 
Due to the low declination of the galaxy, the uv plane 
is unequally sampled with a maximum spacing of 70 m in 
the north-south direction and of 250 m in the east-west. 
The corresponding east-west linear scale at the distance 
of the source, assuming D = 3.4 Mpc (Sandage and Tam- 
mann 1975), is 42 pc. A continuum map was obtained by 
averaging visibilities over a bandwidth chosen to be free 
of line emission. In practice we used spectral channels 
with velocities 200kms~^ above and — 200kms~^ below 
the galaxy's systemic velocity. For the effective 235 MHz 
bandwidth used to obtain continuum data, we derived a 
one a point source sensitivity limit of 0.44mJy/beam. The 
rms noise level in 5 MHz wide channel maps, as determined 
from the dirty maps after subtraction of the continuum, is 
2mJy/beam. This corresponds to an rms brightness tem- 
perature of 17 mK for the synthesized beam size and is 
consistent with a total on-source integration time of 15 
hours and a mean system temperature of 150 K. 

We have estimated the zero-spacing flux lacking in our 
maps. Taking into account that the shortest spacing mea- 
sured by the interferometer is '~20m, we expect to filter 
out scales ^^30". Compared to the 30m data of MP99, 
only ~20% of the single-dish flux in both SiO and H13cO+ 
is missing within the primary beam. Therefore, our data 
offer an unbiased picture of the total SiO and H^^CO"'" 
gas content in the nucleus of NGC 253. 

We assume that the orientation of the NGC 253 disk is 
defined by the angle PAdisfe=51°, and inclination i=78.5°. 
The orientation of the major axis of the stellar bar, along 
PA(,ar=68°, is used on purpose in some figures throughout 
the paper. 

3. The interferometer maps 
3.1. The radio continuum emission 

The 3mm contimium emission contours are plotted in Fig- 
ure la. Despite its compactness, this continuum source, 
with a peak flux of 70mJy/beam, is resolved by our 
beam, with 80% of the flux coming from the central 
source with a deconvolved size of ~8"x4.3", oriented 
along the disk major axis (PA~50°). The integrated flux 
amounts to 0.250±0.008Jy; extrapolating fluxes at 6cm 
and 2cm(Turner and Ho 1983), we estimate that nonther- 
mal and thermal free-free emission account each for one 
third of the observed 3mm continuum flux. The remain- 
ing third would come from emission by dust. Although 



our value for the integrated flux agrees within the errors 
with that measured at 85GHz by Carlstrom et al 1990 
(0.300±0.050Jy), it is slightly lower than that derived by 
P96 at 98GHz (0.320±0.030Jy). 

The better sensitivity of our maps, compared with 
prior observations (Carlstrom et al 1990; P96) reveals 
the presence of two previously undetected sources located 
along the bar major axis at offsets x=-10" and x=10". 
These new sources appear marginally as extensions in the 
1.3mm maps of Kriigel et al 1990. This suggests that they 
correspond mostly to dust emission. However the nature 
of these 3mm new sources is unknown. Comparison of 
our data with measurements at other wavelengths is not 
straightforward because of the missing flux in the images 
at the centimeter range and, on the other hand, owing 
to the lack of resolution of the existent maps at 1.3mm, 
0.8mm and 0.45mm. Before high angular resolution obser- 
vations at short millimeter and submillimeter wavelengths 
are obtained, little can be said on the nature of these new 
sources. 

Contrary to P96, we see no extension of continuum 

emission to the southwest of the nucleus, interpreted as a 
signature of the giant outflow detected in X-rays. Contin- 
uum emission of the dust outflow has only been tentatively 
detected at 0.45mm by Alton et al 1999, but it is absent 
in the 1.3mm and 0.8mm maps. We have extrapolated the 
flux measured at 0.45mm to 3mm and corrected for differ- 
ent resolutions, assuming an emissivity law of . We 
predict a 3mm counterpart of O.lmJy/beam for the dust 
outflow in NGC 253, well beyond the noise level in our 
map (0.44mJy/beam) and a factor of ~200 weaker than 
the 3mm extensions of P96, that reveal as spurious. 

3.2. The SiO and H^^CO^ emission 

Figure 2 shows the velocity-channel maps obtained in the 

(v=0,J=2-l) line of SiO at the centre of NGC 253; figure 
3 shows the same but in the (J=l-0) line of H^'^CO+. Gas 
emission is imevenly distributed in the velocity range [120- 
353kms^^]. This arises from a circumnuclear disk (CND) 
extending over 600pcxl50pc. The disk, seen at high in- 
clination (i=78.5°), is radially resolved and displays an 
overall rotating pattern. 

Although the dominant trend is circular rotation, we 
also flnd emission at highly non-circular velocities in both 
tracers. The emission is patchy in both lines, most no- 
ticeably, the SiO and H^"^CO+ clumps show no one- 
to-one correspondence over the map: as discussed be- 
low (section 7), the T,„5(SiO)/T,„;,(Hi-''CO+) ratio has a 
wide range of variation in the nucleus of NGC 253 (i.e. 
T„6(Si0)/T„b(Hi3C0+) -0.5-8)). 

The CND also shows structure along the apparent mi- 
nor axis. The velocity-integrated intensity maps of SiO 
(Isio) and H^^CO+ {Ihisco+)j shown in Figures Ib-c, 
illustrate the CND substructure. Spatial coordinates are 
arcsecond offsets along the principal axes of the near in- 



4 



S.Garci'a-Burillo et al: Large-scale shocks in the starburst galaxy NGC 253 



Fig. 2. SiO(v=0,J=2-l) velocity-channel maps observed with the IRAM interferometer of Plateau de Bure with a 
resolution (HPBW) of 7.5"x2.6" (PA=180°). Absolute coordinates are 2000.0 (the derived dynamical center is in- 
dicated by the cross at q;j2ooo=00'*47™33.18'', (5j2ooo=-25°17'17.2"). Velocity-channels range from v=361.5kms~^ to 
v=93.8kms~^ by steps of -8.6kms~^. Contour levels are -6, 6, 8, 10, 12, 15, 19 and 23mJybeam~^. Icr-noise level is 
2mJybeam~^. Axes are {a,5). 

Fig. 3. Same as Figure 2 but for the J=l-0 transition of H^"'CO+ with the same levels, velocity spacing and axes 
orientation. 



frared (NIR) stellar bar (Scoville et al 1985). The bar ma- 
jor axis (denoted x) is at PA=68°, with x>0 eastwards 
(the bar minor axis, denoted y. lies along PA=158°, with 
y>0 northwards). Figure lb shows that the SiO emission 
mostly comes from two-nested rings centered at (0,0) and 
roughly oriented parallel to the x axis. A similar morpho- 
logical description applies to H^^CO+ map (Figure Ic), 
although here the inner ring orientation deviates signifi- 
cantly from the bar major axis. 

The inner ring (I), of average diameter D/ ^8"(120pc), 
contains 60-70% of the total flux in the two lines and 
it corresponds to the region of the massive starburst in 
NGC 253 seen in different tracers. The emission of the 
outer ring (II), of average diameter D// ~24"(360pc), 
contains 20-30%of the total flux. The two nested rings 
share as common center the peak of the 3mm contin- 
uum source. Emission in I is detected within a large ve- 
locity range for the two edges of the ring (~120kms~^). 
This spread is symmetrical with respect to the con- 
tinuum source locus (a similar symmetry is seen to 
hold for II). We therefore identify it as the dynamical 
center of NGC 253 (aj2ooo=00''47"33.18" and (5.72ooo=- 
25°17'17.2"). This position coincides, within the errors, 
with the dynamical center determined using optical (Wat- 
son et al 1996), infrared (Sams et al 1994, Boker et al 1998) 
and radio data (Turner and Ho, 1985). 

Although most of the emission arises from the double 
ringed source (I-II), we have also detected gas outside the 
x axis: the CND structure along the minor axis is partly 
resolved. Although part of the emission can be accounted 
by the gas distribution in the CND disk seen in projec- 
tion (NGC 253 is not fully edge-on), at places, the emission 
lies far away from the NIR bar major axis, and most no- 
ticeably, this corresponds to gas clouds showing highly 
non-circular motions and a high relative abundance of 
SiO (T„b(SiO)/T„fe(Hi3CO+)=3-5)). In particular, the 
SiO emission channels of Figure 2 in the range [284- 
335kms~"'^] show receding gas (v-Vsj^s=80kms~-'^; where 
Vsys=235kms^^) in the SE quadrant, where it should be 
approaching. As will be discussed below, the origin of this 
anomalous component is related with the giant outflow de- 
tected in X-ray and optical emission lines. This does not 
have the same origin as I-II. However, there is a H^^CO+ 
clump at (x,y)^(l,8) which shows no SiO counterpart. 
Contrary to the anomalous component, this might corre- 



spond to disk gas seen in projection. This clump has a low 
ratio T„5(SiO)/T^b(Hi3cO+)~0.3, and most important, 
it lies close to the minor axis of the galaxy appearing at 
velocities close to Vsys- 

4. The CND structure 

To study the CND structure we have obtained the 

peak brightness flux distributions in the two lines 
(BP''"''' {x,y)=ma,x[Tmb]v{x,y))- This allows us to filter out 
the lower spatial frequencies enhancing the most intense 
features in the maps. Then we deprojected B^"^"'^ onto the 
galaxy plane, in order to get a sharp picture of the CND, 
viewed face-on. The result is shown in Figures 4a-b. The 
peaks of emission of ring I are rather aligned with the ma- 
jor axis of the disk, tilted 20° clockwise with respect to the 
bar major axis. The tilt is more pronounced for H^^CO+ 
(Figure 4a). On the contrary, the peaks of emission of ring 
II open out, go across the line of the bar major axis and 
delineate a spiral-like arc. The face-on picture of the CND 
of Figure 4a consists of a gaseous spiral- like ridge (II) end- 
ing up in an unresolved ring I viewed edge-on and aligned 
with the major axis of the disk. 

Fig. 4. (a): H^^CO+(1-0) peak-brightness flux contours 
in the CND of NGC 253, deprojected onto the galaxy plane 
((Xg,Yg); where runs parallel to the major axis of the 
disk along PAdisk^51°), displayed from 7 to 19mJy/beam 
by steps of 2mJy/beam. The inner region is a nearly edge- 
on ring encircling the starburst, whereas the outer region 
opens up towards larger PAs, delineating a spiral-like arc. 
(b): SiO(t; = 0, J = 2 — 1) peak-brightness temperature 
contours in the CND of the galaxy, as seen in the plane 
of the sky in (x,y) coordinates as defined in text. Same 
contours as (a). The thick line high-lights the orientation 
of the N and S filaments. 

Interferometer observations by Canzian et al (1988) 
have shown that the bulk of molecular gas in the CND, 
traced by the 1-0 line of ^^CO, is distributed in a bar- 
like source, oriented in the direction of the NIR bar. 
The recent ^'^CO map of Hiittemeister and Aalto (1998) 
confirms this view. In contrast, the distribution of dense 
gas (n(H2)>10^^®cm^'^) is similar to the Bp*^"*' maps of 
SiO and H"CO+. The CN(l-O) and HNC(l-O) maps of 
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Huttemcistcr and Aalto (1998) show a centrally condensed 
inner disk of diameter ^^8", oriented along PA=45°-50°, 
while the outer disk tilts farther towards PA=70°, delin- 
eating a spiral-like arc that goes across the major axis of 
the stellar bar. The same applies for the HCN(l-O) map 
of Paglione et al 1995: the inner disk is better aligned with 
the disk major axis PA=45°, whereas the outer ring tilts 
progressively towards PA=70°. 

There is little doubt that NGC 253 contains a stel- 
lar bar oriented along PA=68° and extended 300" on the 
sky, which sets a lower limit for its corotation radius at 
RcoR=150" (Scoville et al 1985). P96 first explained the 
distribution of the CS(2 1) emission in terms of a bar po- 
tential. However the parameters of the bar were chosen 
arbitrarily: the bar corotation is placed at 40", a factor 
of ^4 too short, compared with the estimated bar semi- 
major axis (~150"). The NIR image of the bar is distorted 
by dust extinction. There are two dust lane ridges offset 
with respect to the bar major axis. Assuming that the 
NW is the near side of NGC 253 (De Vaucouleurs 1958), 
the two offset dust lanes inside the bar represent a trailing 
spiral wave. This pattern is an indication of the existence 
of two Inner Lindblad Resonances (outer=oILR and in- 
ner=iILR; see Athanassoula (1992)). The morphology of 
the gas response in the presence of two ILRs has been 
analysed by several authors (see Buta and Combes 1996 
and references therein). Owing to dissipation by cloud- 
cloud collisions, the gas cannot follow the periodic orbits 
of the stars in a bar (xi=parallel, between corotation and 
oILR or X2: antiparallel, between the oILR and the ilLR). 
The gas response always leads the stellar orbits, giving a 
trailing spiral inside the bar corotation that goes across 
the oILR. It is tempting to identify the spiral-like II as 
the oILR. The gas compression along the spiral ridges 
should be best traced, as observed, by the higher den- 
sity tracers of molecular gas. In contrast, the nuclear spi- 
ral across the oILR is washed out in the ^^CO and ^^CO 
maps tracing the low or moderately dense molecular gas 
(n(H2)=103-4cm-3). 

Moreover, the proposed resonance loci of the bar nat- 
urally explain the double-ringed appearance of NGC 253 
(Figures Ibc), and how the nuclear starburst might have 
been onset. When a barred galaxy has two ILRs, as 
NGC 253, gas tends to accumulate at the ilLR, as a conse- 
quence of secular evolution. The engine behind the radial 
inflow of gas is double. Gravity and viscous torques co- 
operate in bringing the gas closer to the nucleus. Gravity 
torques act on the gaseous spiral expected to form between 
the ilLR and oILR. The spiral is a superposition of lead- 
ing and trailing waves. The leading wave dominates at the 
crossing of the ilLR, whereas the trailing wave develops 
at the crossing of the oILR (Buta and Combes 1996). The 
gravity torques change their sign at the crossing of the 
resonances, and also when we move from the trailing to 
the leading spiral (Combes 1994a, 1994b). The net result 
would be the formation of two pseudo-rings at the ilLR 



and the oILR, with similar time-scales. However, viscous 
torques being always negative, may break this symmetry 
making the gas migrate towards the ilLR in a few dy- 
namical times ('^lO^^^ycars) and therefore progressively 
depopulate the oILR. The oILR may end up as a trailing 
spiral- like pseudo-ring (ring II in NGC 253), whereas the 
leading spiral vanishes into the inner ring (similar to ring 

I)- 

It seems that the morphology of the SiO/H^^CO+ 

maps and how they compare to other molecular gas tracers 
can be explained by gas accumulation at the ILRs of the 
bar in NGC 253. Arnaboldi et al 1995, used the Ha major 
axis kinematics to make an independent estimation of the 
bar resonances loci. Within the epicyclic approximation, 
they identified two ILRs in the nucleus of NGC 253. As- 
suming a lower limit for the bar corotation, Rcor=3.9kpc, 
they derived an oILR at 300pc (20") and the ilLR near 
the center, in agreement with what is proposed in this 
work. Studying the gas kinematics in the SiO/H^'^CO+ 
CND will allow to test the location of the principal reso- 
nances and understand the relation between the starburst 
and the bar. 

5. Kinematics of the CND 

There is a wealth of published work on the dynamics of the 
ISM in the nucleus of NGC 253. The determination of a 
rotation curve (v^ot) for the nucleus of NGC 253 has been 
very controversial because of extinction (Ulrich 1978). The 
high-sensitivity Ha data of Arnaboldi et al 1995, already 
showed a velocity gradient more symmetrical and steeper 
than seen in previous optical surveys. 

Figures 5a-b show the SiO and H^^CO+ position- 
velocity diagrams taken along the kinematical major axis 
of the disk of NGC 253 (x' axis). In view of the high in- 
clination of NGC 253, the terminal velocity method would 
be well suited to derive Vrot from the major axis kinemat- 
ics. However the patchiness of the SiO/H^^CO"'" emission 
makes any fit uncertain. There are strong oscillations of 
the terminal velocities as a flmction of radius, revealing 
an uneven gas distribution and/or the presence of non- 
circular motions. As a zero-order approach we will adopt 
Vrot derived by Canzian et al(1988) from CO. This species 
is the best tracer of molecular gas at a large-scale, from 
densities as low as n(H2)~10^cm~^, and although non- 
circular motions should also affect the observed kinemat- 
ics, the estimated Vrot would be closer to the axisymmetric 
value. 

The ccntroid of SiO/H^'^CO+ emission in the ma- 
jor axis p-v plots follows closely Vrof(CO), but most 
noticeably, the emission shows a large velocity spread 
around (Ai, = 100-120kms~^; see Figures 5a-b). The ob- 
served linewidths could be partly due to the smearing of 
Vrot within our beam. Assuming i=78.5°, we expect to 
get emission from clouds at ydepr=i9-10" (ydepr is the 
deprojected distance from which emission is seen by our 
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Fig. 5. The major axis (along PA=51°) position- velocity diagrams taken at the center of NGC 253 in the SiO(v=0, J=2- 
1) (a, left) and H^"^CO+ (b, right) lines, x' offsets along major axis (PA=51°) are referred to the dynamical center 
derived in this work: aj2ooo=00'*47™33.18^ and (5j2ooo=-25°17'17.2". Vsj,s(LSR) is taken at 235kms"^ Contours are 
-0.004, 0.004, 0.006, 0.008, 0.010, 0.012, 0.015, 0.019 and 0.023Jybcam-i. Our effective resolution along the major axis 
is Ax'x Av=3.2"x8.6kms~^. The rotation curve derived from CO (v,.oi(CO)) is superposed to Fig 5a. The romboid-like 
parallelogram delimits the p-v space occupied by bar-driven orbits inside the ILR domain (see text). 



beam along y'). Along the major axis we get emission from 
x'=±1.6". Vrot (CO) follows a rigid body law with a veloc- 
ity gradient of 6kms^^/", from which we derive an upper 
limit for the total velocity broadening of ~60kms~^ within 
our beam, i.e., a factor of ~2 lower than the observed ve- 
locity spread. 

The existence of a large velocity dispersion in the CND 

gas cannot be excluded (we would need cr^, ^40-50kms~-'^ 
to explain the observed linewidths). However, an explana- 
tion of the SiO/H^^CO+ linewidths in terms of unresolved 
non-circular motions is more plausible as it is supported 
by additional observational evidence. 

Firstly, a part of the emission is detected at velocities 
that cannot be accounted by any circular rotation law. 
These regions in the major-axis p-v plots correspond to 
the (x',v)-quadrants (x'>0,v>V5ys) and (x'<0,v<V5j;s)). 
SiO/H^'^CO"'" gas emission fills unevenly a romboid-like 
parallelogram in Figures 5a-b. The symmetry of the above 
pattern is remarkable and it reminds of the signature 
of bar-driven elliptical orbits developed inside the ILR 
domain (Binney et al 1991; Garcia-Burillo and Guelin 
1995, Fux 1998, Garcfa-Burillo et al 1999). The paral- 
lelogram boundary would correspond to the inner non- 
selfintersecting xi orbit (called the cusped orbit). Inside, 
gas clouds lie along precessing X2 orbits whose envelope 
delineate a spiral across the oILR; these orbits are highly 
elliptical and produce the observed non-circular motions. 
The velocity gradient is reversed for a SiO component in 
the major axis p-v plot between x'=-4" and x'=4". The 
latter was also reported by P96 in the CS p-v plot and and 
by Ananthamaraiah and Goss(1996) using recombination- 
line data. Apparent counter-rotation can also be explained 
by bar orbits (see Figure 8 of Garcia-Burillo and Guelin 
1995). The SiO minor axis position- velocity diagram, dis- 
played in Figure 6a, suggests also that gas distribution in 
the nucleus deviates from axisymmetry. The cut along the 
minor axis is markedly asymmetrical: the centroid of red- 
shifted gas lies l"-2" above the major axis locus, whereas 
blue-shifted gas lies l"-2" below. 

Ananthamaraiah and Goss(1996) reported the exis- 
tence of a complex kinematic subsystem in the inner 150pc 
of NGC 253, using recombination-line data. The nucleus 
would host three gaseous disks; among the anomalous 
components, one exhibits rotation in a plane perpendicu- 
lar to the galactic disk, and the inner disk is apparently 
counter-rotating. These components are introduced to ex- 
plain deviations from circular motions similar to what is 
shown in our data. On the basis of this somewhat ad hoc 



decomposition, the authors find evidence of a past merger 

involving a small mass companion. However, the merger 
hypothesis is unable to explain the undisturbed morphol- 
ogy and kinematics of the outer disk (Olson and Kwan, 
1990). The distorted kinematics of the NGC 253 nucleus, 
including the existence of strong non-circular motions de- 
tected also in the recombination-line data, are more read- 
ily explained in terms of a bar-driven gas flow. 

6. A molecular gas counterpart of the NGC 253 
nuclear outflow 

In the previous section we have argued that the bulk of 

the SiO/H^'^CO"'" emission in the CND fits within the bar- 
driven gas flow scenario. However, we reported in section 
3 the detection of an anomalous component which might 
be far from the CND plane. The off-axis gas, mostly de- 
tected in the SiO line (T„b(SiO)/T„b(H"CO+=3-5), is 
best seen towards the south {down the plane of the CND). 
The southern SiO plume starts at the base of the starburst 
ring I and it reaches a height of (5y=10" (measured from 
the major axis of the bar; see Figure 4b). The gas shows 
also strong deviations from circular motion (A^, ^100- 
150kms~^) and its emission spreads over ~200kms~^ (see 
Figure 6b). There is a much weaker and point-like coun- 
terpart of the southern SiO-plume towards the North. 
The northern gas shows also strong noncircular motions 
(A„ '~-150kms~^) but it is here characterized by a narrow 
profile (^17kms~^). SiO spectra towards these positions 
are displayed in Figure 7. The Northern and Southern 
plumes (hereafter called N and S) can be connected by a 
line that goes across the dynamical center; the N-S line 
has a PA=131° measured from North (see Figure 4b). 
This orientation places the anomalous component close 
to the galaxy minor axis (PA=141°) where no contribu- 
tion from rotation curve in the radial velocity gradient 
is expected. Moreover it is hard to explain non-circular 
motions of A^, ^100-150kms~^ due to bar forcing in the 
plane of the galaxy. 

There is compelling evidence for the existence of a gi- 
ant outfiow of gas in the nucleus of NGC 253. Demoulin 
and Burbidge(1970) and Ulrich (1978) were the first to re- 
port the presence of an ionized gas outflow, based on the 
analysis of gas kinematics derived from Ha spectroscopy. 
Soft X-ray emission from the giant outflow was later dis- 
covered by Fabbiano and Trinchieri (1984). The X-ray 
nebula is more extended to the southern part of the galaxy 
minor axis and was interpreted as thermal emission of gas 
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Fig. 6. a SiO strip along the minor axis of the galaxy (PA=141°) with same contours as Figure 5. The effective 
resolution is Ay'x Av=4"x8.6kms^^. We also display in b the SiO position-velocity diagram parallel to the minor 
axis, but taken at x=4", showing the presence of gas emission at highly non-circular motions (inside the high-lighted 
region). 

Fig. 7. Wc compare the SiO(v=0,J=2-l) (thick line) and H^''CO+ (thin line) spectra at several positions: 
(x,y)=(0",0")(left bottom) in the dynamical center, (x,y) = (3.2",3.2")(left middle) in region I, (x,y)=(- 
17.6",1.6")(left top) in region II, (x,y)=(3",-8") (right top) in the southern outflow and (x,y) = (-5",8")(right 
bottom) in the northern outflow; the spectra illustrate the large variations of the Trnb{S'^0)/TmbOi^^CO~^) ratio 
in the CND. Y axis scale is in Jy/beam. 



heated by fast shocks, driven by a bipolar wind coming 
from the starburst. Further support for this scenario is 
provided by the distribution of ionized gas in the outflow, 
which adopts the form of fllaments that border the X-ray 
nebula (McCarthy et al 1987). Furthermore, Schulz and 
Wegner (1992) detected line-splitting in Nil, SII and HII, 
consistent with the ionized gas lying across the conical 
working surface of the X-ray nebula. The gas entrainment 
by the hot wind may cause shock-heating and subsequent 
emission. The observed line-splitting, together with the 
inferred expansion velocities (~390kms~^), are the main 
arguments supporting the shock-heating scenario. 

Very recently, Alton et al 1999 have reported the ten- 
tative detection of submillimeter emission at 450/u coming 
from one fllament associated with the S outflow. The lo- 
cation of the S filament agrees reasonably well with the 
one detected in SiO/H^^CO"*" (see above and Figure 4b). 
The dust mass contained in the filaments, which partly 
depends on the dust temperature, cannot be easily es- 
timated from the detection at just one wavelength. No 
counterpart of the dust outflow is seen neither at 1mm 
nor at 3mm, as expected from the flux detected at 450/x 
(see discussion of section 3.1; see also Kriigel et al, 1990). 
However wc can try to estimate the dust cohimn density 
towards the S filament from the flux reported by Alton et 
al 1999, corrected to our beam. We assume a dust emis- 
sivity spectral index of (3=2, the emissivity law is taken 
from Chini et al (1997), and adopt a dust temperature 
in the range Tfi„s(=13K-37K (similar to the model that 
fits the M82 dust outfiow: Alton et al 1999). We cal- 
culate a range of Ndust ~0.6-4x lO'^M© towards the S 
filament, within our beam. If Mgas/Md„st=100, this re- 
sults implies N(H2)~0.6-4x lO'^Mo/b cam. As discussed 
in section 7, we can estimate N(H2) in the S filament 
where H^^cq-f detected, using N(Hi3C0+) derived 
from a Large Velocity Gradient (LVG) transfer model 
and assuming X(H^'^CO+)=10~^° (see discussion of sec- 
tion 7 for details). We calculate N(H2)'~3.5x lO^Mg/beam 
(4.7x lO^^cm^^), roughly in agreement with the value es- 
timated from the dust emissivity. Therefore no contradic- 
tion exists between the failed detection of the S outflow in 
the continuum and the successful detection of the outflow 
in the two lines of SiO and Hi3C0+. 



If the molecular CND is locally disrupted and it is en- 
trained by the hot wind out of the plane, we can expect 
that the geometries for the working surfaces of molecu- 
lar and ionized gas should be alike. The parameters of 
Schulz and Wegner's model are an expansion velocity of 
390kms^^ for the conical outflow coming out of the plane, 
and a half-opening angle of 25°. Projection effects along 
the line of sight would make appear two velocities at 
each position of the outflow (N or S). The emission from 
the S(N) fllament should appear either at Vrad=-Wkms^^ 
(120kms-i) or at Vrad=300kms-i (440kms-i). The veloc- 
ities observed in the SiO fines for the N (vrad=100kms~^) 
and S filaments (300kms~^) confirm the predictions of 
Schulz and Wegner's model, originally intended to fit the 
observed line-splitting in the optical lines. Turner (1985) 
was the first to report the existence of a molecular gas 
ejection event in the NGC 253 nucleus. However, the lo- 
cation of the OH plume hardly fits the morphology and 
kinematics of the ionized or molecular outflow: the OH 
gas is ejected from the nucleus towards the NE quadrant 
(along PA=6°) and it reaches ~lkpc height (see Fig 11 of 
Turner(1985) for details). 

We conclude that the distribution and kinematics of 
the off-axis component detected in our maps are easily 
accounted for if we suppose that the emission comes from 
molecular gas in an outflow leaving the nucleus, out of the 
galaxy plane. Furthermore, an estimate of the molecular 
gas mass involved in the outflow seems to be consistent 
with the dust emissivity of the filaments. 

7. The T^6(SiO)/T„5(Hi3cO+) and 
N(SiO)/N(Hi3cO+) ratios 

The line temperature ratio 

RTmb=T„fe(SiO)/T„b(Hi3CO+) shows large varia- 
tions (0.5-8) in the nucleus of NGC 253; these variations 
however show a systematic trend. For the sake of sim- 
plicity, we first examine the ratio of velocity-integrated 
emission Rj=I(SiO)/I(H^^CO+) averaged across those 
regions where R/ shows a characteristic value. We have 
chosen as cases of study the dynamical center, regions 
I-II and the outflow components (N-S). Table 1 shows 
the mean values of R/ and its range of variation; Figure 
7 shows representative spectra for the regions studied. 
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Wc sec a clear trend in the ratios, which increase as 
one moves away from the vicinity of the starburst. The 
average ratios are R/ ~1.0, for the center, Rj ~1.1 for 
the starburst ring I, R/ '--^3 for the outer region II and 
R/ ~2 for the outflow. We will analyse the range of 
variation of Rrmfi for each region defined above, and 
derive the column density ratios N(SiO)/N(H^^CO+) via 
a Large Velocity Gradient (LVG) scheme. Our objective 
is to study the range of physical parameters and chemical 
abundances of SiO clouds in the CND. 

LVG solutions depend on three parameters: the den- 
sity (n(H2)), the kinetic temperature (Tk) and the abun- 
dance of the species. The kinetic temperature of molec- 
ular gas in the center of NGC 253 has been derived by 
several authors using multitransition studies in ^^CO and 
its isotopes (13C0 and C^^O), CS and HON; all indi- 
cate an average value of T^ >50K. We have adopted 
Tfe=50K as representative for the SiO clouds and taken 
n(H2) from the fit to the ratios R32=I(SiO(3-2))/I(SiO(2- 
1)) and R54=I(SiO(5-4))/I(SiO(3-2)) derived by MP99, 
using single-dish 30m spectra. MP99 found little evidence 
of a density decrease for the SiO clouds, at least within 
the inner r~20" of the nucleus. They derived an aver- 
age n(H2)=5xlO^ assuming a common filling factor for 
the three lines. Wc will take this estimation as represen- 
tative for SiO and H^'^CO+ in our model. Note that al- 
though there might be a hotter component in the nucleus 
(Tfe ~100K; see Mauersberger et al 1990), we must stress 
that n(H2) inferred from the SiO line ratios depend weakly 
on the adopted kinetic temperature if >50K. 

Figure 8 displays the range of LVG solutions that 
fit R32 and R54 within the errors. The assumed density 
n(H2)=5xlO^ fits satisfactorily R32 and R54, for a large 
interval of N(SiO)/At,. Trying to delimit the LVG solution 
along the N/A„ axis is equivalent to choosing a beam fill- 
ing factor {rjfiii) for the clouds. Whereas it is justified to 
assume that the two lines share a common jy/iH, its exact 
value in an external galaxy is unknown a priori. It is worth 
discussing to what extent our conclusions might depend 
on riftu . 

We have run LVG models for n(H2)=5x lO^cm"'^ and 
Tfc=50K, fitting RT„fe=T„(,(Si0)/T„6(Hi3C0+) in the 
different regions of the CND. Table 1 fully explores the de- 
pendence of N(SiO)/N(H"CO+) on rjfui; the filling factor 
is varied within its whole range (0.005-1). The main con- 
clusion, independent of rjfm, is that the SiO abundance 
is significantly enhanced in the outer CND (region II), 
where we derive N(SiO)/N(Hi3CO+)=5-15, and contrary 
to what it might be expected, it is lower at the center and 
at the starburst ring I, where N(SiO)/N(H"CO+)=l-2. 
The SiO abundance rises in the outflow, where we infer 
N(Si0)/N(Hi3C0+)=3-5. 

A density fall-off with radius has been reported 
by Wall et al 1991 (and conflrmed by Harri- 
son et al 1999). They concluded that n(H2) might 
reach ~5xl0^cm~^ for r>20", not far from re- 



Fig. 8. This figure represents LVG solutions fit- 
ting the R54=I(SiO(5-4))/I(SiO(3-2)) and R32=I(SiO(3- 
2))/I(SiO(2-l)) line ratios (continuous lines) observed in 
the SiO clouds of the nucleus of NGC 253, assuming 
T=50K. We consider errors of 20% and 30% in ratios R54 
and R32 (pointed and dashed lines delimit the region of 
allowed solutions within the errors tolerance). An aver- 
age density of n(H2)=5xl0^cm~'^ is a good compromise, 
indicating that the SiO emitting gas is dense. Trying to 
delimit further the LVG solution along the N/A„ axis is 
equivalent to choose a beam filling factor {rjfui) for the 
clouds. 

gion II. If we take n(H2)~5xlO*cm~^, we find that 

N(SiO)/N(Hi-''CO+)=20-100 for II. Opacity in the 2 1 
line of SiO is close to ~1 for this solution, whereas r ~0.02 
for H^^CO"'". To fit the ratio of line temperatures we 
need imperatively to rise N(SiO)/N(H^''CO+) in region 
II. Consequently, conclusions on the SiO abundance en- 
hancement of region II remain unchanged, if not clearly 
favoured, in the lower density scenario. 

From the ratio N(SiO)/N(Hi3CO+), we can infer the 
absolute abundance of SiO (X(SiO)=N(SiO)/N(H2)), if 
we adopt a canonical value for H^^CO+ (<X(H^^CO+)>). 
It is expected that SiO will be more sensitive than 
H^^CO"'" to any variation of the chemical environment 
in NGC 253. This hypothesis is successfully confronted to 
the measurements of abundances of H^^CO"'' obtained in 
a large variety of molecular clouds in our Galaxy where 
values close to 10~^ are obtained in all cases (from which 
we derive X(H"CO+)-'10~io ff we adopt i^c/i^c^l/gO). 
Though there is observational and theoretical evidence 
that X(Hi2co+) may be reduced to 10"^° in the hot 
cores of GMCs, SiO is lower by a similar factor in these 
sources (see Blake et al 1987); moreover, very compact 
sources are not expected to dominate the emissivity of 
H^^CO"'"(1-0) at the scale we are observing the nucleus 
of NGC 253. Contrary to X(Hi3C0+), X(SiO) is mea- 
sured to vary by several orders of magnitude in the 
Galaxy, depending on the type and the location of molec- 
ular clouds in the disk (see discussion). We can estimate 
<X(Hi3C0+)> in the nucleus of NGC 253. From our data 
we obtain a column density of <N(II^^C0+)>=3x 10^^, 
averaged within the CND. Similarly, we can calculate 
<N(Il2)> from the ^^CO(1~0) map of Canzian et al 
(1988), averaged on the same area. If we take the con- 
version factor x=N(H2)/Ico=3xl0^°cm2Kkm-is, it fol- 
lows <N(H2)>~3x 10^^, which implies global abimdanccs 
of <X(H"CO+)>~10-io in NGC 253, quite similar to the 
value inferred for our Galaxy. 

Column 5 of Table 1 illustrates the variation of X(SiO) 
in the CND of NGC 253. We derive X(SiO)~l-2xlO-^° 
for the center and the starburst ring I, X(SiO)~3- 
15x10-1° for region II, and X(SiO)-3-5xlO-i° for the 
outfiow. The global abundance of SiO in the CND would 
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Table 1. Tabic 1 fully explores the dependence of N(SiO)/N(H-'^^CO"'") and X(SiO) on rjfiu ; the filling factor is varied 
within its whole range (0.005-1) and the results are displayed for the different regions of the nucleus where distinct 
values of R/ are measured. See text for details. 



be <X(SiO)>~1.5xlO namely, an order of magnitude 
above the typical value for a PDR (Janssen et al 1995, 
Walmsley et al 1999). Most notably, X(SiO) reaches the 
highest abundance in the outer region II, where the chem- 
ical processing of the gas by the nuclear starburst, mostly 
adscribed to region I, should be minor. The SiO abun- 
dance is also significantly enhanced in the outflow, cer- 
tainly powered by the starburst, but not directly on it. 
This result emphasizes that SiO emission is not always 
directly related to the dense gas near sites of recent star 
formation. Therefore, different mechanisms must be ex- 
plored to account for the different chemical processing of 
molecular gas within the nucleus. 

8. Discussion 

We have established the existence of three distinct regions 
in the nuclear region of NGC 253, based on morphologi- 
cal, kinematical and chemical criteria. Region I dominates 
the global SiO emission and it is associated with the nu- 
clear starburst taking place across ~150pc. It presents 
a ring-like pattern, interpreted as the gas response near 
the ilLR of the bar. The SiO abimdance in the starburst 
(X(SiO)~l-2xlO~^°) is significantly larger than that es- 
timated for a typical PDR (X(SiO-PDR)~10-ii). Region 
II extending up to the map edges (r^300pc(20")); dis- 
plays a spiral-like morphology also present in other high- 
density tracers. The observed distribution and kinemat- 
ics of molecular gas arc interpreted as the signature of 
the bar oILR. The orbit crowding across the oILR coin- 
cides with the detection of strong non-circular motions 
(Av^50 lOOkms^^) and, most notably, with an enhance- 
ment in the SiO abundance (X(SiO)~a few 10^^). The 
third region constitutes the molecular counterpart of the 
nuclear gas outflow observed in X-ray and optical lines. 
This was recently identified as a dust chimney. Two fila- 
ments (N and S) come out of the plane of the NGC 253 
nucleus from r±60pc, near a working surface where the 
gas is entrained by the outflow. Again, we have found here 
a link between the existence of high-velocities (Av~100- 
150kms-i) and an SiO enrichment (X(SiO)>3- 5x 10"^°). 

Two processes are thought to dominate the chemistry 
of the molecular material associated with a starburst: the 
photodissociation of the gas by intense UV field produced 
by OB associations, and the existence of strong shock 
waves. Shocks can be generated locally during mass loss 
episodes in young stars (often identified as bipolar out- 
fiows). Alternatively, shocks may take place on a larger 
scale, either produced by supernova explosions or related 
to spiral/barred density waves. Contrary to X(H^^CO"'"), 
X(SiO) is measured to vary by several orders of mag- 



nitude in the Galaxy, depending on the type and the 
location of molecular clouds in the disk: star forming 
clouds (X(SiO)~10~^-10^^), galactic center clouds appar- 
ently not forming stars (X(SiO)~10~^) or quiescent clouds 
(X(SiO)<10-i2-)_ 

It is well established that shocks occurring in bipo- 
lar outflows can enhance the abundance of SiO in the gas 
phase to reach, in the most extreme cases, ~10~^ (Martin- 
Pintado et al 1992). Although the starburst region I might 
be described as a giant PDR (Carral et al 1994), the rel- 
atively large abundances of SiO measured and the high 
densities of SiO clouds indicate that most of the emis- 
sion arises in bipolar outflows powered by young massive 
stars, similarly to what is observed in the galactic com- 
plexes W51 and W49. Moreover, the particularly large ra- 
tio of I([SiII])/I([OI])~l measured in the inner region of 
NGC 253 indicates an enhancement of Si in gas phase. 
This was interpreted by Carral et al 1994 as evidence of 
grain destruction in the starburst region, most likely by 
shocks related to massive star formation. 

The case for shock processing is even stronger in region 
II and in the outflow where X(SiO)~(50-100)xX(SiO- 
PDR). However, the origin of the large abundance of SiO 
in region II is less clear. SiO emission could arise partly 
in bipolar outflows powered by young massive stars like 
in region I. The major drawback to this explanation is 
that region II is located away from the nuclear starburst. 
Therefore one should not expect an enhancement of SiO 
when the chemical processing caused by star formation 
has declined. Alternatively, the SiO enhancement could 
be produced by large-scale shocks in the molecular gas, 
caused by the crossing of clouds orbits near the oILR 
of the bar. The only evidence for chemical enrichment of 
gas attributed to large-scale shocks is found in the Galac- 
tic Center where values of ~10~^ have been measured in 
SiO clouds (Martm-Pintado et al 1997, Hiittemeister et al 
1998). However there is no theoretical consensus whether 
the shocks can be caused by the action of density waves, 
especially in molecular clouds. The hierarchy of fragmen- 
tation inhibits shocks in molecular clouds; the mean-free 
paths for cloud-cloud collisions are of the same order of 
magnitude as the thickness of the disturbing potential well 
(Casoli and Combes 1982, Combes and Gcrin 1985). 

If the starburst scenario gives rise to SiO emission in 
bipolar outflows of young massive stars, the differences in 
the abundance of SiO in regions I and II might be ex- 
plained if the starbursts are in a markedly different stage 
of evolution. The older starburst would occur in region I, 
where a large fraction of the young massive stars are al- 
ready on the main sequence. The latter explains why the 
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classical tracers of the starburst phenomenon are mostly 
restricted to I. In contrast, region II would contain a 
younger and less-evolved starburst, where most of the 
young massive stars would be still in the mass loss phase, 
characteristic of the pre-main sequence stage, where ener- 
getic bipolar outflows occur. 

Finally, the interpretation of the morphology and kine- 
matics of the molecular gas filaments give additional sup- 
port to the outflow scenario. The enhancement of SiO in 
the gas phase observed in this region can be explained 
by the chemical processing of grains by the strong high- 
velocity shocks (v>200kms~^) associated to the outflow. 
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